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ABSTRACT 

We have studied the radio and far-infrared (FIR) emission from 114 galaxies in the 7 nearest clusters 
(< 100 Mpc) with prominent X-ray emission to investigate the impact of the cluster environment on the 
star formation and AGN activity in the member galaxies. The X-ray selection criterion is adopted to 
focus on the most massive and dynamically relaxed clusters. A large majority of cluster galaxies show 
an excess in radio emission over that predicted from the radio-FIR correlation, the fraction of sources 
with radio excess increases toward cluster cores, and the radial gradient in the FIR/radio flux ratio is 
a result of radio enhancement. Of the radio-excess sources, 70% are early-type galaxies and the same 
fraction host an AGN. The galaxy density drops by a factor of 10 from the composite cluster center out to 
1.5 Mpc, yet galaxies show no change in FIR properties over this region, and show no indication of mass 
segregation. We have examined in detail the physical mechanisms that might impact the FIR and radio 
emission of cluster galaxies. While collisional heating of dust may be important for galaxies in cluster 
centers, it appears to have a negligible effect on the observed FIR emission for our sample galaxies. The 
correlations between radio and FIR luminosity and radius could be explained by magnetic compression 
from thermal ICM pressure. We also find that simple delayed harassment cannot fully account for the 
observed radio, FIR, and mid-IR properties of cluster galaxies. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (A262, AWM7, A426, 
A1060, A1367, Virgo, A1656) — galaxies: elliptical and lenticular, cD — galaxies: 
interactions — infrared: galaxies — radio continuum: galaxies 



1. INTRODUCTION 

The radio and far-infrared (FIR) properties of local 
galaxies with a wide range in bolometric luminosity, from 
the faintest dwarf galaxies to ultraluminous infrared galax- 
ies (ULIRGs), have been well documented. The over- 
whelming majority of these galaxies obey the linear corre- 
lation between radio and FIR luminosity which holds over 
5 orders of magnitude in luminosity (Price & Duric 1992). 
The mechanisms for radio and FIR emission are thought 
to be well understood and both have been shown to be ac- 
curate tracers of recent star formation (see Yun, Reddy, & 
Condon 2001-hereafter YRC-and references therein). The 
small percentage of luminous galaxies which deviate from 
the FIR-radio relation in the local universe are thought 
to be powered by AGN (for those showing excess radio 
emission over that expected from the radio-FIR correla- 
tion and excess mid-IR emission) or are in a short-lived 
dust enshrouded AGN/starburst phase (for those show- 
ing excess FIR emission; YRC). Deviations from the cor- 
relation may also arise from environmental changes such 
as the presence of an intracluster medium (ICM). Studies 
have found notably lower FIR-radio flux ratios for cluster 
galaxies than expected from the FIR-radio relation in the 
field (see Miller & Owen 2001a), indicating some modi- 
fication of the radio and/or FIR emission mechanisms in 
cluster galaxies. 

Mechanisms that can only affect the observed FIR emis- 
sion include collisional heating of dust grains by ICM 
electrons and ram pressure stripping of ISM. In the lat- 



ter case, observations indicate that such stripping only 
strongly affects the outer HI disk of the galaxy, leaving the 
dense molecular gas and dust component intact (Kenney 
& Young 1986, Cayette et al. 1994). A consequence of the 
lower star formation rates observed in cluster galaxies (e.g, 
Gomez ct al. 2002) is a corresponding reduction in both 
FIR and radio emission. Enhancement of radio continuum 
and FIR emission may result from tidally induced star for- 
mation. Mechanisms affecting the observed radio emission 
include an increase in the magnetic field strength via gas 
compression as galaxies move through the ICM (Gavazzi 
& Boselli 1999; Rengarajan, Karnik, & Iyengar 1997; An- 
derson & Owen 1995; Scodeggio & Gavazzi 1993), simple 
thermal compression of the magnetic field (e.g., Miller & 
Owen 2001a), and/or elevated AGN activity arising from 
frequent galaxy-galaxy collisions or galaxy harassment in 
the cluster environment (e.g., Moss & Whittle 1993; Lake, 
Katz, & Moore 1998). Studying radio and FIR emission 
should therefore be useful in examining these mechanisms. 

Miller & Owen (2001a) conducted a study of the radio 
and FIR properties of a sample of 329 local Abell cluster 
galaxies. Optically-identified AGNs in their sample show 
greater radio-to-FIR flux ratios with more scatter than 
star forming galaxies and galaxies in cluster cores show 
enhanced radio emission. While unable to rule out ram 
pressure enhancement, they find that thermal compression 
of the magnetic field alone could produce the observed 
radio luminosity segregation. 

Prior investigations into the effects of the ICM on galaxy 
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properties have primarily selected clusters based on their 
optical properties that may not accurately reflect the state 
of the intracluster medium. For example, the Miller & 
Owen (2001b) cluster sample has 329 spectroscopically 
confirmed cluster galaxies, half of which reside in poor 
clusters where the environmental effects are not as pro- 
nounced as in rich clusters. We have adopted a different 
approach in selecting clusters based on the presence of X- 
ray emission. The most common explanation for this emis- 
sion in clusters, first proposed by Felten et al. (1966), is 
that X-rays result from thermal bremsstrahlung from hot 
(« 10 8 K), low density (~ 10~ 3 atoms cm~ 3 ) intraclus- 
ter gas. The inclusion of X-ray luminous clusters in our 
sample (1) ensures that only massive, dynamically-relaxed 
clusters are chosen, and (2) is immune to the biases in- 
herent in optically-selected surveys, such as sensitivity to 
morphological type, AGN fraction, and ages and masses 
of galaxies. 

X-ray cluster and radio/FIR galaxy selection allows us 
to compile a homogeneous sample of cluster galaxies in 
order to address the key questions of this study: (1) does 
the rich cluster environment help or hinder star formation 
activity in cluster galaxies, and (2) how does the pres- 
ence of an ICM affect the frequency of AGN sources? To 
aid in addressing these questions, we also examine the fol- 
lowing: convolution due to projection effects (§3), clus- 
tering of early-type sources as predicted by the density- 
morphology relation (Dressier 1980; §6.1), and FIR emis- 
sivity as a function of mass (§6.2). All of these effects can 
influence the observed radio and FIR properties of cluster 
galaxies and must be carefully disentangled. To determine 
whether rich clusters preferentially affect low or high lu- 
minosity galaxies, we also examine the faint and bright 
end behavior of the cluster luminosity function (LF) and 
compare with that found in the field. 

The organization of this paper is as follows. In §2, we 
present the criteria used to select clusters and individual 
galaxies. The effects of projection are quantified in §3. 
The radio-FIR relation and flux density ratio for cluster 
galaxies are discussed in §4. The effects of the ICM on the 
galaxy radio and IR luminosity functions are discussed in 
§5. The radio and FIR properties of the sample arc in- 
terpreted in the context of the density-morphology rela- 
tion, mass segregation, and various emission mechanisms 
in §6, including electron-dust collisional heating and ram- 
pressure and thermal magnetospheric compression. We 
also examine the mid-IR properties of the cluster sample 
to address the frequency of AGN sources in §6. 

2. SAMPLE SELECTION 

2.1. Cluster Selection 

The clusters used in this study were culled from the 
flux-limited sample of 55 bright X-ray clusters compiled 
by Edge et al. (1990) using the following selection criteria. 
First, we chose only those clusters with 5 > —40°, corre- 
sponding to the survey limit of the NRAO VLA Sky Sur- 
vey (NVSS) from which radio data were obtained. Clus- 
ters near the galactic plane < 10°) were also excluded. 
Finally, a redshift cutoff of z < 0.025 was applied (100 
Mpc for n o = 75 km s _1 Mpc~ : — this value of H G is as- 
sumed throughout the paper), to ensure accurate NVSS 
and IRAS detections given the sensitivity and confusion 



limits of the two surveys. As an added advantage, we only 
select nearby, well-studied clusters, so that a case-by-case 
study of individual galaxies can be done. This selection 
process yielded 7 clusters: A262, AWM7, A426 (Perseus), 
A1060 (Hydra), A1367, Virgo, and A1656 (Coma). The 
ACO (Abell, Corwin, & Olowin 1989) cluster richness 
among these 7 clusters ranges from -1 to 2, but all have 
virial masses of ~ 10 14 M to ~ 10 15 M Q (David, Jones, 
& Forman 1996; Koranyi & Geller 2000; Ettori, Fabian, 
& White 1998; Tamura et al. 2000; Donnelly et al. 1998; 
Collcss & Dunn 1996; Schindler, Binggeli, & Bohringer 
1999) indicative of a well-defined potential. The central 
position, assumed redshift, and 2 — 10 kcV X-ray luminos- 
ity of each cluster are included in the data summarized in 
Table 1. 

2.2. Galaxy Selection 

Individual galaxies in the 7 aforementioned clusters 
were selected using the following procedure. We used the 
NASA/IPAC Extragalactic Database (NED) to search for 
all cataloged galaxies with available redshift within 1.5 
Mpc projected distance from the cluster center. Sources 
with radial velocities > 2000 km s _1 than the systemic 
cluster redshift were excluded from the sample. NED se- 
lection is subject to the availability of data, and to inves- 
tigate this bias we applied the selection criteria above to 
the Miller & Owen (2001b) sample for clusters common 
between their sample and ours: A262, A426, A1367, and 
Coma. The same was done with the AWM7, A1060, and 
Virgo homogenous samples of Koranyi & Geller (2000), 
Richter (1989), and Binggeli, Sandage, & Tammann (1985), 
respectively. All of these samples reduce to subsets of our 
sample after imposing our selection criteria. There are 3, 
2, and 5 galaxies in A262, A1367, and Coma, respectively, 
which are not included in the homogeneous samples above, 
and may be drawn from other surveys having potentially 
different selection biases. However, we note that these 
sources have radio and infrared (IR) luminosities compa- 
rable to the rest of the sample, and arc morphologically 
diverse (half are early- and half are late- type galaxies), so 
their inclusion in the sample should not significantly affect 
our conclusions. 

A constant velocity cutoff of 2000 km s _1 preferentially 
excludes sources with high peculiar velocity, particularly 
for A1060 and Coma, both having line-of-sight velocity 
dispersions of a ~ 1000 km s _1 . Adopting a 3cr cutoff 
velocity for A1060 and Coma adds only 1 source with de- 
tected radio and FIR flux to the sample, the starburst 
galaxy NGC4858. Our velocity cutoff quantitatively only 
affects the normalization of the cluster luminosity func- 
tion (§5), since there is evidence that all types of galaxies 
fill luminosity space uniformly (i.e., no mass segregation; 
§6.2). 

The FIR and radio data were obtained from the In- 
frared Astronomical Satellite (IRAS; Neugebauer et al. 
1984) and the NRAO VLA Sky survey (NVSS; Condon 
et al. 1998), respectively. The IRAS data include 12/zm, 
25^m, 60/um, and 100/zm flux measurements with a lcr er- 
ror ellipse of 15". The total luminosity between 40/im and 
120/im, dominated by young massive stars, can be esti- 
mated from the 60/xm and 100/iin flux. The NVSS radio 
survey at 1.4 GHz (20 cm) with the D and DnC configu- 
rations of the Very Large Array, has an angular resolution 
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Table 1 
Cluster Properties 



Name 


a 

(J2000.0) 


<5 

(J2000.0) 


< cz > a 
(km s- 1 ) 


(km s- 1 ) 


Class c 


Nftjre 

(Obs) 


(Obs) 


^early 




2-10 kcV L x d 
(10«h- 2 75 ergs- 1 ) 


A262 


01 52 50.4 


36 08 46 


4887 


588 





6 


7 


5 


8 


1.22 


AWM7 


02 54 32.2 


41 35 10 


5168 


680 


-1 


1 


4 


2 


3 


5.20 


A426 


03 18 36.4 


41 30 54 


5366 


1324 


2 


1 


5 


1 


5 


48.9 


A1060 


10 36 51.3 


-27 31 35 


3777 


647 


1 


5 


16 


4 


17 


1.09 


A1367 


11 44 29.5 


19 50 21 


6595 


879 


2 


11 


9 


8 


12 


3.08 


Virgo 


12 26 32.1 


12 43 24 


1075 


632 


1 


3 


20 


2 


21 


0.52 


Coma 


12 59 48.7 


27 58 50 


6925 


1008 


2 


6 


20 


11 


15 


33.3 



a Systemic cluster velocity from NED 

b Velocity dispersion (la) from Strubble & Rood (1999), Koranyi et al. (1998), and Fadda et al. (1996) 
c Richness class from Abell, Corwin, & Olowin (1989) 
d 2 - 10 kev luminosity from Edge et al. (1990) 



of 45" with a la sensitivity of 0.5 mJy beam 1 at 20 cm. 
The NVSS is sensitive to galaxies with star formation rates 
;> 1.4 M Q yr" 1 within a distance of 100 Mpc (YRC; Con- 
don 1992). 

In the case of the two nearest clusters, A1060 and Virgo, 
we searched the NVSS database for radio sources within 10 
kpc projected distance (the typical size of a disk galaxy) 
from the optical position of each of the galaxies. For the 
other 5 clusters, we chose a search radius of 30" based 
upon the NVSS and IRAS error ellipses. The probabil- 
ity of finding an unrelated NVSS source within the search 
radius is < 0.1%. In some cases, confusion between a 
foreground galaxy and a background radio source (e.g., 
UGC1416), or between interacting companions (NGC4567 
and NGC4568), was resolved using a more detailed in- 
spection. For nearby galaxies with large angular extent in 
Virgo, the NVSS underestimates the flux due to extended 
disk emission and in some cases produces multiple radio 
matches for a single optical position. In these cases, cor- 
rections to the radio flux were obtained from the multiple- 
snapshot study of Condon (1987). 

IRAS SCANPI coadditions of each galaxy were obtained 
and the source is included in our sample if the IRAS source 
centroid is within 30" of the nominal radio position from 
the NVSS database and has S/N > 3 in the 60^m and 
100/im bands. The IRAS satellite has a raw resolution 
of ~ 4' and for extended sources greater than this, the 
SCANPI software can miss extended flux. There are 24 
such sources in Virgo. IRAS fluxes for 14 of these extended 
sources were obtained from the IRAS Bright Galaxy imag- 
ing compilation of Soifer et al. (1989). Accurate fluxes 
could not be found in the literature for 10 of the extended 
objects in Virgo, but all of these sources are excluded from 
the sample after adopting a luminosity cutoff (see §2.3). To 
summarize, we selected galaxies from NED with rcdshifts 
and they are included in the sample if they have secure 
NVSS and IRAS 60/xm and 100,uni detections. This selec- 
tion process yields 182 galaxies in the composite cluster. 

The 1.4 GHz radio and 60^m luminosities were calcu- 



lated using the following relations: 

l0gil.4GHz(W Hz" 1 ) = 20.08 + 21og£> + logSi.4GHz (1) 

logL 6OMm (L ) = 6.014 + 21og£> + log5 60 ^ m (2) 

where D is the distance in Mpc, and Si.4GHz and Smiim 
are flux densities in units of Jy. The 60^m luminosity is 
related to the FIR luminosity, Lfir, by 

Lfir(Lq) = ( 1 + *°° Mm ) L 60t im (3) 

\ /.OOD60/jm/ 

(see Helou et al. 1988). The typical uncertainty in the 
luminosities is ~ 15% (oiogL ~ 0.065), based on NVSS 
and IRAS absolute flux calibration errors. Intrinsic scatter 
due to variations in FIR and radio emission dominates the 
uncertainties in the flux scaling (YRC). 

2.3. Sample Limits and Divisions 

We have imposed a luminosity limit of logL 60Mm = 8.92 
(the 3a IRAS detection limit for the most distant clus- 
ter in our sample, Coma, at 92 Mpc). Figure 1 shows 
the distances for the galaxies in the 7 clusters as a func- 
tion of 60^m luminosity. Over 37% of the sample is ex- 
cluded by adopting our luminosity- limit, including the 10 
extended Virgo galaxies for which accurate IRAS fluxes 
could not be found. The final luminosity-limited catalog 
contains 114 sources. The sample is small compared to 
optically-selected cluster samples, but is drawn from a ho- 
mogeneous selection of clusters based on X-ray emission 
and is luminosity-complete. The catalog of 114 galaxies 
is summarized in Table 2, along with the galaxies' blue 
magnitudes and RC3 morphological codes from NED. 

A key motivation for this study is to understand the 
effect of ICM density on the star forming activity in mem- 
ber galaxies. Using radius as a proxy for local density, 
each cluster has been divided into a Core and Ring region. 
Core galaxies are defined as having a projected cluster- 
centric distance less than 0.5 Mpc (see §3 for a justifica- 
tion of this radius), and Ring galaxies are defined to have a 
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Table 2 

Cluster Galaxies in the Luminosity-Limited Sample 







a 


5 


m 


< cz > 


log Lx.AGHz 


log L60 M m 






Cluster 


Name 


(2000.0) 


(2000.0) 


(mag) 


(km s~ 1 ) 


(W Hz" 1 ) 


(L ) 


q 


Notes 


A 9R9 
,rt_ZOZ 




1 AA zll 


Qr: 17 ir; i 
OO ±1 ±0.± 


±o.oo 


A~\ ^1 
4101 


91 fiS 
Zl .Oo 


y. ID 


9 ^7 






UGC1308 


1 50 51.34 


36 16 34.3 


13.77 


5171 


22.01 


9.00 


1.29 


E 




UGC1319 


1 51 28.99 


36 3 56.7 


14.50 


5310 


21.95 


9.63 


2.03 


Irr 




NGC703 


1 52 39.30 


36 10 16.0 


14.27 


5580 


21.63 


9.09 


1.91 


so- 




UGC1347 


1 52 45.98 


36 37 8.1 


13.49 


5543 


21.79 


9.82 


2.36 


SAB(rs)c 




NGC708 


1 52 46.45 


36 9 8.3 


13.70 


4855 


22.53 


9.10 


0.98 


cD;E;Sy2 




ivrPP7i n 

IN t LU 


1 ^9 9^ 
1 DZ D4.ZO 




1 A. 97 
14. Z ( 


9^ 
01ZD 


99 9^ 
ZZ.ZO 


Q £9 

y .oz 


1 .04 


OCQ 




KUG0151+356 


1 54 11.72 


35 54 57.0 


15.70 


4387 


21.11 


8.99 


2.36 


Spiral 




UGC1385 


1 54 53.81 


36 55 4.3 


13.90 


5621 


21.98 


10.43 


2.64 


(R)SB0/a;Sbrst 




UGC1416 


1 56 46.18 


36 53 3.3 


14.72 


5470 


21.11 


8.97 


2.17 


Sb 




NGC753 


1 57 41.99 


35 54 58.8 


12.97 


4903 


22.23 


10.20 


2.31 


SAB(rs)bc 




NGC759 


1 57 50.43 


36 20 37.2 


13.84 


4667 


21.89 


9.60 


2.00 


E 




UGC1460 


1 59 4.61 


36 15 25.2 


14.64 


4874 


21.24 


9.50 


2.53 


Sa 


AWM7 


IC259 


2 49 45.80 


41 3 7.7 


15.00 


6108 


21.50 


9.74 


2.52 


SB0? 




NGC1106 


2 50 40.49 


41 40 17.8 


13.30 


4337 


22.88 


9.78 


1.10 


SA0+;Sy2 




UGC2350 


2 52 40.27 


41 23 45.8 


14.60 


3964 


21.22 


9.38 


2.51 


Sb 




NGC1122 


2 52 50.76 


42 12 18.1 


12.90 


3599 


21.80 


10.03 


2.53 


SABb 




CGCG539-121 


2 54 16.65 


42 43 34.1 


14.70 


6331 


21.91 


9.82 


2.18 


Spiral 


A426 


UGC2608 


3 15 1.48 


42 2 8.8 


13.70 


6998 


22.76 


10.63 


2.06 


(R')SB(s)b;Sy2 




UGC2617 


3 16 1.32 


40 53 14.7 


13.80 


4697 


21.33 


9.53 


2.69 


SAB(s)d 




IC310 


3 16 42.95 


41 19 28.7 


13.89 


5660 


23.01 


9.50 


0.86 


SA(r) 




UGC2654 


3 18 43.38 


42 17 59.9 


14.20 


5793 


21.87 


9.94 


2.38 


S? 




NGC1275 


3 19 48.16 


41 30 42.1 


12.64 


5264 


25.15 


10.52 


-0.47 


cD;pec;NLRG;Sy2 




UGC2696 


3 22 5.79 


42 10 15.3 


15.50 


5453 


21.25 


9.09 


2.13 


Scd 


A1060 


IRAS10288-2824 


10 31 13.55 


-28 39 51.8 


15.33 


3690 


20.98 


9.20 


2.43 


SABa 




ESO436-G034 


10 32 44.35 


-28 36 40.2 


14.38 


3624 


21.21 


9.32 


2.50 


Sb 




ESO501-G017 


10 34 24.09 


-26 29 27.8 


14.97 


4429 


21.17 


9.46 


2.47 


S0(7) 




NGC3285B 


10 34 37.22 


-27 39 12.5 


13.86 


2952 


21.09 


9.08 


2.39 


SB(rs)b 




ESO436-IG042 


10 34 38.74 


-28 35 0.3 


14.40 


3190 


21.77 


9.87 


2.28 


E4 




ESO436-G046 


10 34 50.35 


-28 35 4.7 


13.42 


3436 


20.92 


9.60 


2.97 


SB(rs)b 




ESO437-G002 


10 34 59.31 


-28 4 44.2 


14.63 


2311 


21.30 


9.23 


2.28 


Sc;pec 




ESO437-IG003 


10 35 7.89 


-27 59 27.2 




2391 


21.22 


8.95 


2.15 








in ^ 9^7^ 

1U OO io. t o 


OO IO AQ Q 


1 °. Q°, 
ij.yj 


oou*± 


91 0^ 


9 21 


9 ^ 






ESO437-G015 


10 36 57.88 


-28 10 31.2 


13.53 


2743 


20.99 


9.09 


2.32 


SB(s);sp 




NGC3312 


10 37 2.55 


-27 33 52.7 


12.68 


2886 


22.17 


9.30 


1.56 


Sab(r)p;AGN 




ESO501-G045 


10 37 12.64 


-26 40 22.5 


15.11 


4578 


20.92 


9.17 


2.52 


S0(3)/a(r) 




NGC3314A 


10 37 12.64 


-27 41 3.1 


14.40 


2872 


21.91 


9.60 


1.99 


Sab:sp 




ESO501-G053 


10 37 37.91 


-26 16 37.0 


14.80 


3814 


21.06 


9.33 


2.56 


Sa 




ESO501-G059 


10 37 49.46 


-27 7 17.4 


14.97 


2437 


21.66 


9.36 


1.99 


Sc 




ESO501-G065 


10 38 33.53 


-27 44 12.4 


13.71 


4425 


21.86 


9.55 


2.00 


SB(s)d:pec 




ESO437-G025 


10 38 40.10 


-28 34 4.9 


14.49 


3458 


20.98 


9.19 


2.57 


S(B)b 




IRAS10369-2827 


10 39 16.78 


-28 42 52.1 


15.69 


2811 


21.12 


9.17 


2.24 


SO/a 




ESO501-G068 


10 39 17.47 


-26 50 23.7 


14.25 


3090 


20.93 


9.15 


2.57 


S(B)bc(rs) 




NGC3336 


10 40 17.15 


-27 46 46.4 


13.00 


4000 


21.72 


9.59 


2.22 


SA(rs)c 




ESO501-G075 


10 40 58.75 


-27 5 8.8 


13.48 


4750 


21.18 


9.33 


2.56 


SA(s)c 


A1367 


CGCG097-068 


11 42 24.49 


20 7 11.8 


14.55 


5974 


22.05 


10.15 


2.35 


Sbc 




CGCG097-073 


11 42 56.22 


19 58 10.3 


15.50 


7275 


22.30 


9.37 


1.29 


SAcd:pec 




UGC6680 


11 43 2.27 


19 38 36.1 


15.04 


6985 


21.79 


8.94 


1.71 


Sb 




CGCG097-079 


11 43 12.22 


20 26.8 


15.70 


7000 


21.67 


9.38 


2.03 


Irr 




UGC6697 


11 43 48.45 


19 58 14.3 


14.08 


6725 


22.70 


10.10 


1.66 


Im 




CGCG097-092 


11 43 58.49 


20 11 16.3 


15.50 


6373 


21.43 


9.52 


2.33 


Sbc 




NGC3840 


11 43 59.28 


20 4 32.5 


14.54 


7368 


21.60 


9.79 


2.45 


Sa 




IN uL'Oo4Z 


I 1 AA 1 Q1 

II 44 l.Ol 


iy DO DO. f 


1Z. < O 


OolO 


99 nfi 
zz.uo 


Q AA 
y .44 


i. ( y 


Hi 




IRAS11419+2022 


11 44 32.05 


20 6 22.0 


15.40 


7230 


21.46 


9.94 


2.70 


SO 




UGC6719 


11 44 47.47 


20 7 25.5 


14.44 


6571 


21.42 


9.30 


2.15 


Sab 




NGC3860B 


11 44 46.80 


19 46 13.3 


15.30 


8293 


21.73 


9.81 


2.45 


Irr 




NGC3860 


11 44 48.86 


19 47 47.1 


14.22 
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1.98 


Sa 




NGC3859 


11 44 51.20 


19 27 21.5 


14.76 


5468 


21.97 


9.90 


2.24 


Irr 




CGCG097-125 


11 44 54.97 


19 46 24.7 


15.74 


8271 


21.35 


9.69 


2.71 


SOa 




NGC3861 


11 45 4.33 


19 58 21.0 


13.47 


5085 


21.91 


9.51 


2.07 


(R')SAB(r)b 




NGC3862 


11 45 5.23 


19 36 37.8 


13.67 


6511 


24.51 
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9.72 
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11 46 12.23 
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9.37 


1.65 


SA(r)0/a;LINER;Syl 
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Table 2 Continued. 



Cluster Galaxies in the Luminosity-Limited Sample 
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projected clustercentric distance between 0.5 Mpc and 1.5 
Mpc. We chose the maximum limit of 1.5 Mpc since this 
is a median value for the virial radii of most clusters, typ- 
ically between 1 Mpc and 3 Mpc. Some fraction of Core 
galaxies are Ring galaxies appearing in projection and the 
effect of this contamination is addressed in §3 and §4.3. 

Morphological dependencies are investigated by divid- 
ing the sample into early-type (E and SO) and late-type 
subsamples. Of the 182 galaxies, 54 (30%) are early-type 
galaxies, higher than the detection rate of early-type FIR 
galaxies in the field (20%; see Cotton & Condon 1998), 
and significantly higher than the fraction of ellipticals in 
the local universe (~ 15%). This clustering of early-type 
galaxies reflects the density-morphology relation (Dressier 
1980; §6.1). It may also signify that cluster ellipticals are 
brighter at FIR wavelengths than those found in the field 
(§6.1). The number of Core, Ring, early-, and late-type 
galaxies in the luminosity-limited sample for each cluster 
are summarized in Table 1. 
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Fig. 1. — Distribution of sources as a function of 60/im luminosity. 
Each horizontal band is composed of all detected galaxies within a 
particular cluster. The dashed and solid lines denote the flux and 
luminosity limit (logZ/60pm = 8.92), respectively, adopted for the 
sample. Galaxies below the luminosity limit comprise 37% of the 
flux-limited sample, and most of these are in Virgo. 

3. DEPROJECTION AND DENSITY DISTRIBUTION OF 
CLUSTER SOURCES 

Accurate information on the nature of the Core galax- 
ies can be obtained by estimating the amount of source 
contamination from the Ring region due to projection ef- 
fects. A geometrical recursive technique employing sim- 
ple discretization of the Abel integral is used to deproject 
the sources within spherical shells surrounding the core, 
starting with the outermost shell. The parameters that 
characterize the density distribution are the central den- 
sity (p ), core radius (r c ), and log-slope (/3; see equations 
4 and 5 below), and of these, p can be the most system- 
atically uncertain: the outer shell has the lowest number 
count, and this statistical error can become significant af- 
ter propagating to the inner shells, resulting in a range of 



normalizations. In contrast, (3 and r c are more resilient to 
such initial error as they parameterize a "bulk" behavior 
of the sources. Minimizing the initial error in p a requires 
accurate background decontamination, in which case ve- 
locity selection may be used to determine the most likely 
cluster members, as is done in this study. The clustcr-to- 
cluster variation in density distribution is dominated by 
statistical noise, but there is a clear trend in that of the 
composite cluster. Table 3 lists the radii of the spherical 
shells, observed counts within those regions, computed de- 
projected counts, and estimated density, which increases 
with decreasing radius. The "radius" of a bin is defined as 
the mean of the inner and outer radii of that bin. 

In relaxed, regular clusters, the number density of galax- 
ies and the gas density are correlated since the crossing 
time of a cluster galaxy is an order of magnitude shorter 
than a Hubble time if the cluster is in hydrostatic equilib- 
rium. The hydrostatic model for these clusters provides 
a consistent fit for both the X-ray and number density 
distributions (Bahcall & Lubin 1994). The galaxy surface 
density profile can be represented as 

/ r 2 \ (-3/2)0+1/2 

n = n (l + -j (4) 

according to the /3-model (Cavaliere & Fusco-Femiano 
1976). The deprojected density that should fit both the 
number and X-ray gas distribution is then given as: 

/ r 2N ("3/2)/3 

P = Po[l + ^J • (5) 



55 I 1 1 1 1 1 1 1 1 1 1 r 




I i i i I i i i I i i i 

0.5 1 1.5 

r (Mpc) 

Fig. 2. — Observed galaxy surface density (open circles) with best- 
fitting surface density fi profile of equation (4) (solid curve). Our 
best-fit values are /3 ~ 1.05, r c ~ 0.64 Mpc, and n a ~ 49 Mpc -2 , 
with x 2 = 2.09. The density drops by a factor of ~ 10 from Mpc 
to 1.5 Mpc. 

The accuracy of the deprojection method can be deter- 
mined by fitting the galaxy surface density with equation 
(4) (Figure 2). While the 95% confidence limits for n Q , r c , 
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Table 3 
Deprojection Statistics 





Radial Boundaries 






Number Density 


Region 


(Mpc) 


Observed Counts 


Deprojected Counts 


(Mpc- 3 ) 


(Core) 


0.0 < r < 0.5 


33 


15.92 


30.41 


1 


0.5 < r < 0.7 


20 


17.85 


19.55 


2 


0.7 < r < 0.9 


17 


10.16 


6.28 


3 


0.9 < r < 1.1 


18 


19.49 


7.73 


4 


1.1 < r < 1.3 


16 


22.47 


6.19 


5 


1.3 < r < 1.5 


10 


28.11 


5.70 


Total 




114 


114.00 





and (3 are a factor of two times the best-fit values, the mean 
values of n and r c are consistent with those found for typi- 
cal clusters (n G = 49 Mpc" 2 , r c = 0.64 Mpc, and /3 = 1.05 
with x 2 = 2.09). Comp arison with the independent fit 
for the deprojected galaxy density (not shown) yields pa- 
rameters that agree with those in Figure 2. These results 
indicate that selection based on radio and FIR emission 
does not bias the galaxy distribution away from that seen 
among clusters in general. 

4. FIR-RADIO RELATION 
4.1. Radio-FIR Correlation 

Figure 3 shows the 20 cm radio continuum luminosity 
plotted against the 60/xm luminosity, i.e., the radio-FIR 
correlation, for the flux-limited sample. The band between 
the two dotted lines indicates the best-fit linear correlation 
between radio and FIR luminosity and the scatter of 0.26 
dex for the 1809 galaxies of the IRAS 2 Jy YRC field sam- 
ple. Most (98%) of the IRAS selected field galaxies, includ- 
ing early- types (E and SO), are normal star forming galax- 
ies following the linear radio-FIR correlation (YRC). Com- 
parison with these field galaxies shows that ~ 37% of the 
cluster galaxies deviate significantly (> 0.26 dex) from the 
radio-FIR correlation defined by the IRAS-selected sam- 
ple. For comparison, < 2% of field galaxies deviate > 0.26 
dex from the mean relation, mostly from systematic devi- 
ations at the high and low luminosity ends. Twenty-seven 
percent of the cluster galaxies above our flux limit exhibit 
a higher radio luminosity than expected from the linear 
field correlation (> la above the field correlation): 50% 
of these galaxies are in the observed Core region and 62% 
are early-type galaxies. 

None of the sources have Lqq^ > 10 11 L Q , in contrast 
to the field where ~ 11% are luminous in the infrared 
(YRC). The 2 Jy selection at 60^m is identical to the op- 
tical selection of late- type galaxies in the field (YRC). We 
expect the same selection effect for our cluster sample. 
Rarity of IR luminous galaxies has been noted in other 
similarly-sized samples of cluster galaxies (e.g, Bicay & 
Giovanclli 1987). From YRC, the space density of IR lu- 
minous galaxies in the field is about 10~ 5 Mpc -3 , and 
there is only a 0.1% chance of finding such a source within 
the search volume (~ 100 Mpc 3 ) in the field. The high 
galaxy density in the cluster environment should increase 
the expected fraction by an order of magnitude to near 
unity. Their non-detection in these clusters by itself is not 



a strong indication of bias against IR luminous phenom- 
ena in the cluster environment. For comparison, the Abcll 
cluster study of Miller & Owen (2001b) finds only 4 of 329 
cluster sources with secure IRAS 60/im and 100/im detec- 
tions having Lfir > 10 11 L©, and only 1 of these 4 lies 
> 1.5 Mpc from the cluster center. 

A separate linear least-squares fit to the radio and FIR 
luminosity for our 182 sample galaxies indicates a shal- 
lower than unity log-slope of 0.76, but with a dispersion 
of 0.74, emphasizing the much greater scatter among the 
cluster galaxies as compared with those from the field (Fig- 
ure 3) . The scatter is slightly less if we only consider galax- 
ies in the luminosity-limited sample. The shallower slope 
is not statistically significant given the large scatter among 
cluster galaxies. 

Note the apparent bias of sources below the best-fit cor- 
relation between logLgo^m = 8.00 and our luminosity cut- 
off of logL60^m = 8.92. Approximately 90% of the galax- 
ies below the luminosity cutoff are in Virgo. The excess 
of radio-luminous galaxies at logL 6 o M m = 7.5 L Q consists 
of mostly elliptical and/or Scyfert (9/16) galaxies, none of 
which are included in the Virgo radio-FIR study of Niklas, 
Klein, & Wielebinski (1995). One source (IC3099) in this 
region is one of the 10 for which we could not find accu- 
rate IRAS fluxes: the radio bias of this source may be due 
to missing extended FIR emission. Niklas et al. (1995) 
identify a few Virgo early-type galaxies in the Core with 
excess 2.8 and 6.3 cm radio emission compared with their 
FIR emission. This is also the case for early-type galaxies 
further away from the cluster center (Figure 3). Excluding 
these radio-luminous objects, the majority of the remain- 
ing galaxies are situated below the radio-FIR correlation 
for the YRC field sample, and a formal fit to the data yields 
a greater than unity slope, consistent with the Niklas et 
al. (1995) study and with the nonlinearity seen at lower 
luminosities in the radio-FIR correlation of the YRC field 
sample (Fig. 5a of YRC; see also Bell 2002 for a more 
recent review) . Possible explanations for this nonlinearity 
include dust heating by low-mass stars or cosmic-ray dif- 
fusion. Adopting a luminosity cutoff excludes this nonlin- 
ear regime and selects the region in luminosity-luminosity 
space where both the FIR and radio luminosity are directly 
proportional to the massive star formation rate, allowing 
for an unbiased method of segregating star forming galax- 
ies and AGNs (see §6.3). 
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Fig. 3. — Plot of 20 cm continuum radio emission (eq. [1]) versus 
IRAS 6Ctytm luminosity (eq. [2]) for the flux-limited sample. The 
sloping solid line denotes the best-fit linear correlation for the YRC 
field sample and the two dotted lines indicates its la dispersion. The 
vertical line delineates our adopted luminosity-limit. Large circles 
identify early-type and Core galaxies in Figure 3a and Figure 3b, 
respectively, which are biased above the field radio-FIR correlation. 
The small points refer to late-type and Ring galaxies in Figure 3a 
and Figure 3b, respectively. Note the significant number of sources 
above the correlation compared with the few that fall below. 
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Fig. 4. — Logarithmic FIR-radio ratio, q parameter, plotted for 
each galaxy as a function of IRAS 60/im luminosity, with separate 
emphasis on morphological and radial dependence in Figure 4a and 
Figure 4b, respectively. Filled circles and open squares represent 
early- and late-type galaxies, respectively, in Figure 4a. Filled grey 
triangles and circles represent early- and late-type Ring galaxies, 
respectively, and open squares denote Core galaxies in Figure 4b. 
The solid horizontal line denotes the average (q) = 2.34 found for 
YRC field galaxies, dotted lines delineate IR- and radio-excess re- 
gions, and the solid vertical line denotes our luminosity-limit of 
logL 6 o Mm = 8.92. 
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4.2. FIR-Radio Flux Density Ratio 

To investigate the deviation from the linear radio-FIR 
correlation, we compute the q parameter (Condon, Ander- 
son, & Helou 1991), 



q = log 



FIR 



3.75 x 10 12 W mT 2 



log 



Si 



4GHz 



W m^Hz- 1 



(6) 



where S± 4GHz is the observed 1.4 GHz flux density in units 
of W m~ 2 Hz' 1 and 

FIR = 1.26 x Kr 14 (2.58S 60Aim + S w ^ m ) W m" 2 , (7) 

where and Sioo^m are IRAS 60 /xm and 100 jum band 

flux densities in Jy (see Helou et al. 1988). The q parame- 
ter is both independent of distance and starburst strength 
(YRC; Lisenfeld, Volk, & Xu 1996) and so is a good quanti- 
tative measure for comparing galaxies independent of their 
star formation rates. The q values for the cluster galaxies 
are listed in Table 2 and are plotted in Figure 4. The lines 
appearing at q = 3.04 and q = 1.64 delineate the point at 
which the galaxies have 5 times larger IR and radio flux, 
respectively, than what would be expected from the radio- 
FIR correlation. Sources with q > 3.04 and q < 1.64 are 
referred to as IR- and radio-excess sources, respectively 
(YRC). 

The mean q for the flux-limited sample is 2.07±0.74, sta- 
tistically consistent with the mean value of q = 2.34 ±0.01 
for the field. However, the low q value and large scat- 
ter as compared with the field is due in part to AGN 
activity, which results in a systematic difference in the q 
distribution. The mean q for the luminosity-limited sam- 
ple is similar, (q) = 2.10, with smaller dispersion, despite 
the exclusion of low radio luminosity sources among the 
low IR luminosity population of Virgo discussed above. 
This indicates that the radio-luminous sources below our 
luminosity cutoff (e.g., elliptical and Seyfert galaxies in 
Virgo) balance the effect of the galaxies appearing below 
the radio-FIR correlation. About 1/2 of the Abell cluster 
galaxies of Miller & Owen (2001a) populate poor clusters 
and are projected up to 3 Mpc from the cluster center, yet 
have a lower mean (q) = 1.95 (a = 0.77). This is likely due 
to their selection based strictly on radio emission and the 
presence of a large number of optically-identified AGNs in 
their sample. 

4.2.1. cD Galaxies 

There are a total of 13 radio-excess sources in our sam- 
ple, 4 of which are cD galaxies (NGC708, NGC1275, 
NGC4839, NGC4874). Radio and X-ray observations of 
these cD galaxies imply jet-like activity in their nuclei 
(Conselice, Gallagher, & Wyse 2001; Holtzman et al. 1992; 
Vermeulen, Readhead, & Backer 1994; Rhee, Burns, & 
Kowalski 1994; Neumann et al. 2001; Colics & Dunn 
1996; Feretti & Giovannini 1987). The FIR luminosity 
may result from star formation, but in some cases, such as 
NGC1275, there is also a notable contribution from colli- 
sional heating of dust (e.g., Irwin, Stil, & Bridges 2001 — 
see §6.4.1 for further discussion). 

4.2.2. Remaining Radio-Excess Sources 

Of the remaining 9 radio-excess sources, 5 show radio, 
X-ray, and/or emission lines consistent with the presence 



of an AGN. (Bonatto et al. 1996; Owen, Ledlow, & Keel 
1996; Sambruna, Eracleous, & Mushotzky 1999; Veron- 
Cetty & Veron 1996; Ji et al. 2000; Baum et al. 1997; 
de Koff et al. 2000; Bai & Lee 2001; Zirbel & Baum 
1995). CGCG097-073 is a spiral galaxy in A1367 that 
shows strong and extended Ha emission, indicating vigor- 
ous star formation activity (Moss, Whittle, & Pesce 1998). 
KUG1258+287 is an HI rich spiral galaxy in Coma with 
relatively low optical surface brightness in the Digital Sky 
Survey image and may be interacting with its close com- 
panion FOCA195 (Bravo- Alfaro et al. 2000). NGC4927 
is another early-type Coma galaxy. NGC3884 is close to 
the boundary for a radio-excess galaxy, with q = 1.65. It 
is a Seyfert galaxy, with an extended LINER spectrum 
(Durret et al. 1993). 

Thirteen of 114 galaxies (11%) show radio-excess, 11 
times higher than the detection rate in the YRC field sam- 
ple. Of the 13 radio-excess sources, 9 are early-type galax- 
ies. Of the 15 and 18 early-type galaxies in the Core and 
Ring region (in the luminosity-limited sample), 40% and 
17%, respectively, are radio-excess sources. As discussed 
above, at least 5 of the non-cD radio-excess sources host 
an AGN. The average q for both late- and early-type ob- 
served Core galaxies is 1.74, a FIR-radio flux density ratio 
half that of the cluster average. The average q in the Ring 
is 2.24, still a factor of 1.3 times higher than the mean 
q for the YRC field sample. These results suggest some 
variation in q with radial distance. 

4.3. q-Radius Correlation 

The distribution of galaxies as a function of their q val- 
ues is shown in Figure 5 (binned every 0.25 dex). Objects 
with negative q values are included in the first bin. The 
distribution for Ring galaxies peaks near the mean q value 
for the YRC field sample, whereas the peak for the Core 
galaxies lies almost a factor of 2.3 lower. Both distribu- 
tions drop off quickly for high q values, and taper off more 
slowly for lower q values, indicating a general bias toward 
lower q values. 

The deprojection analysis (Table 3) indicates that 17 of 
the observed 33 Core members could be Ring galaxies seen 
in projection. Assuming these contaminating sources have 
the same q distribution as the Ring sources, then 95% of 
the observed Core galaxies with q > 2.00 are Ring mem- 
bers, and will be removed after deprojection and "added" 
to the Ring sample. If this statistical result reflects a real 
change in the population, then the average q value of the 
Core galaxies drops from 1.74 to 1.18, a factor of 14 de- 
crease from that of field galaxies (YRC). 

Of the Ring galaxies, 83% have q > 2.00. Similarly, 
91% of Core galaxies have q < 2.00. If we assume that 
all galaxies with q < 2.00 are Core galaxies and all with 
q > 2.00 are Ring galaxies, then "deprojecting" accord- 
ing to q value indicates that 14 and 95 galaxies are in the 
Core and Ring regions, respectively. For comparison, de- 
projection based on the number counts yields 17 and 99 
Core and Ring sources, respectively. If the above assump- 
tions are made, then the q value can be used to roughly 
discriminate between Core and Ring galaxies, and would 
be useful in situations where the density profile of cluster 
galaxies cannot be determined accurately. A more mean- 
ingful segregator may be observed with a larger sample 
of data and should also vary with core radius and cluster 
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Fig. 5. — Distribution of q values for the luminosity-limited sam- 
ple of cluster galaxies. Sources with q < are included in the first 
bin. The hashed and unhashed portions of the distribution are for 
observed Core and Ring galaxies, respectively. The solid vertical 
line delineates the average q value for the YRC field sample. De- 
projection analysis suggests that most of those Core galaxies with 
q > 2.00 arc Ring galaxies seen in projection (see text). 



richness. A detailed investigation is necessary to establish 
the strength of, and variation in, the q segregator. 

According to Figure 4, 7 of 33 (21%) observed Core 
galaxies have radio-excesses, and 6 of these are early-type 
galaxies. Similarly, 6 of 81 (7%) Ring galaxies have radio- 
excesses, and 2 of these are early-type galaxies. The in- 
crease in frequency of radio-excess sources toward the cen- 
ter of clusters contributes to a radius-dependent q distribu- 
tion. There is significant scatter in the correlation between 
q and r, and the best fit relation is 



q = (0.43 ± 0.12) x r + (1.84 ± 0.10), 



(8) 



where the 4 cD galaxies in our sample have been excluded 
for this analysis — all of which are classified as radio-excess 
sources based on our definitions — since these are not typ- 
ical of cluster galaxies. 

If all radio-excess sources are excluded, the correlation 
between q and r is shallower with a scatter of 0.28 dex, 
identical to the scatter in q for YRC field galaxies. Con- 
sequently, the scatter docs not reflect statistical noise, but 
is a systematic result of the radio-excess sources in our 
sample, most of which are AGNs. If the increased ra- 
dio luminosity of elliptical galaxies stems from an increase 
in number density, then the correlation of equation (8) 
reflects the density-morphology relation as the early-type 
galaxies have a density distribution that peaks in the Core 
region. 

Figure 6 shows the radio and FIR luminosity (computed 
from eq. [3]) for the luminosity-limited sample (exclud- 
ing cD galaxies) as a function of radius. Formal fits to 
the correlations give logLi.4GHz °c (—0.45 ± 0.14) x r and 
logieo^m oc (-0.01±0.09)xr, giving q cx (-0.01 + 0.45) x 
r = (0.44 ± 0.23) x r, and indicates the general bias of ra- 
dio luminous sources toward the Core region. We also 



note that the fraction of galaxies with logLi.4GHz < 21.00 
increases from 0% in the Core to 18% in the Ring. FIR 
luminosity is constant with radius. The lack of FIR lu- 
minosity segregation has also been noted by Rengarajan, 
Karnik, & Iyengar (1997) for spiral cluster galaxies. We 
find the same to be true for early-type galaxies. Hence, q 
values are lower in the Core than the Ring region due pri- 
marily to an increase in the radio luminosity, as opposed 
to a deficit in the FIR luminosity. Possible variations in 
the FIR luminosity-radius correlation between clusters are 
investigated in §6.4.2. 
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Fig. 6. — logLgrjfim and logLi.4GHz luminosity as a function of 
radius for the luminosity-limited sample in Figure 6a and Figure 6b, 
respectively. The solid lines indicate the best-fit correlations. The 
dashed curve indicates the expected radio luminosity as a function 
of radius for a field spiral galaxy undergoing thermal compression by 
the ICM (see §6.5.2). Solid circles are early-type galaxies and open 
squares are late-type galaxies. The correlation for the IRAS 60/im 
luminosity is flat, whereas the radio luminosity shows some variation 
with radius. The correlation between q and r can be completely 
accounted for by the correlation between logL^QUz an d r - 
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5. RADIO AND IR CLUSTER LUMINOSITY FUNCTIONS 

Comparison of the cluster galaxy luminosity distribu- 
tion with an appropriately chosen field sample should draw 
out distinctions between the cluster and field populations, 
and we can thereby investigate whether the rich ICM pref- 
erentially affects the faint or bright ends of the cluster 
galaxy population. We do this analysis by computing the 
cluster radio and IR luminosity functions (LFs) via the 
1/V method (Schmidt 1968; Felten 1976)— i.e., number of 
sources per luminosity interval divided by the total vol- 
ume equal to the number of clusters falling below the lu- 
minosity interval times the volume of one cluster. For 
simplicity, we have assumed the volume of a single cluster 
to be that of a sphere of radius 1.5 Mpc. In practice, the 
sampled volume (and hence normalization of the LF) may 
be difficult to constrain due to a number of factors (e.g., 
projection effects, the adoption of a 2000 km sec -1 veloc- 
ity cutoff, etc.). The deprojection results of §3 suggest 
that the assumption of a spherical volume is not unrea- 
sonable. Furthermore, our adopted velocity cutoff should 
only affect the normalization of the LF since all types of 
galaxies (e.g., Core, Ring, early-, and late-type) are found 
to uniformly populate the range of luminosities considered 
here. We are most interested in comparing the shapes of 
the cluster and field LFs and not the exact normalization. 
Finally, a joint analysis of the radio and IR LFs should 
mitigate some of the uncertainties introduced by sample 
selection and volume determination. 

The LF is constructed by first computing the number of 
galaxies per luminosity interval, taken to be A log L 60fnn = 
0.25 for logLgo^m between 7.50 and 11.00. The cluster LF 
is compared to the best-fit Schechter (1976) LF for the low- 
luminosity YRC field population, normalized to the value 
of the LF in the luminosity bin containing the greatest 
number of sources (representative of L* galaxies). 

5.1. IRAS 60/im Cluster Luminosity Function 

The cluster LF is well-described by the shape of the low- 
luminosity field LF (Figure 7; YRC) for log L mtlla > 9.0, 
after renormalizing the field LF by a factor of ~ 180. The 
break in the LF at logL 60tim <~ 8.5 is a result of the seven 
clusters not being uniformly distributed in redshift space, 
where the sampled volume increased by a factor of 2. For 
comparison, the LF computed for Virgo cluster galaxies 
is also shown in Figure 7. The paucity of IR and radio 
sources in the Virgo core (§6.4.2) indicates that the sam- 
pled volume for Virgo may be smaller than the one taken 
here, yet correcting for this does not fully account for the 
disparity between the cluster and field LFs at low lumi- 
nosity. The bias of the cluster LF below the field LF may 
therefore be a result of incompleteness resulting from non- 
detections of low luminosity Virgo sources. Previous stud- 
ies indicate a ~ —1.3 for optically-selected Virgo galaxies 
(e.g., Ferguson & Sandage 1988). Assuming that the FIR 
emission from faint blue galaxies is correlated with their 
optical luminosity, the faint-end slope for the Virgo clus- 
ter LF in this study is inconsistent with the Ferguson & 
Sandage (1988) result. The divergent behavior of the low- 
luminosity end of the IR LF is discussed below (§5.3). 

5.2. 1.4 GHz Cluster Radio Luminosity Function 



Figure 8 shows the cluster radio LF, constructed in the 
same way as the IR LF. The cluster radio LF is scat- 
tered around the field LF. Unlike the situation with the 
IR LF, the low-luminosity cluster radio sources appear to 
have a faint-end slope similar to that of the field. This 
can be easily explained by the increased contribution from 
radio- luminous Virgo sources with logLgo^m between 7 
and 8 (e.g., Figure 3). There is an excess of cluster galax- 
ies at the high luminosity end which is more pronounced 
than that seen in Figure 8, as there are 3 sources with 
logLi.4GHz > 24.5 that are not shown in the figure. The 
sources showing excess radio emission are seen more clearly 
in Figure 4, where the scatter of cluster galaxies is not 
symmetric about the scatter among YRC field galaxies. A 
similar excess is found in the radio LF for YRC galaxies, 
and after removing all sources within 1.5 Mpc of the cen- 
ters of known clusters within the limits of the YRC field 
sample, the field galaxies have a LF indistinguishable from 
the UGC sample LF of the Condon (1989) study. The ex- 
cess observed in Figure 8, as well as the excesses observed 
for the IRAS 2 Jy sample (see Figure 8b in YRC) and 
the radio LF of UGC sample galaxies (Condon, Cotton, 
& Broderick 2002) all suggest that radio-luminous AGNs 
play a more significant role in the cluster environment than 
they do in the field. 

If radio-excess is an indicator of AGN activity (see §6.3), 
then the excess observed in Figure 8 is predominantly due 
to the early-type radio galaxies, which account for ~ 70% 
of the radio-excess objects (cf. §4.2). To quantify this 
further, the radio LF for late-type galaxies only is com- 
puted and shown in Figure 8. While early- type galaxies 
contribute to the luminosity density for the entire lumi- 
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Fig. 7. — The IRAS 60/im LF for the 182 sources in the composite 
cluster is shown by the open circles. The solid curve denotes the 
best-fit low luminosity field Schechter LF for the 60/im data of the 
YRC field sample, normalized to the value of the cluster LF for the 
bin containing the greatest number of sources. The systematic bias 
of the cluster LF below that of the best-fit field correlation suggests 
incompleteness at the low luminosity end. For comparison, we have 
also shown the IR LF for Virgo cluster sources only, denoted by the 
filled triangles, along with its best-fit luminosity function indicated 
by the dashed curve. 
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nosity range considered here, the difference between the 
LF for the late-type galaxies and the LF for the aggre- 
gate sample increases for log.Li.4GHz > 22.0, suggesting an 
increasing contribution of radio emission from early-type 
galaxies at high luminosities. The 3 galaxies for which 
logii.4GHz > 24.5 are early-type Core galaxies. Despite 
the uncertainty involved in normalizing the cluster LF, no 
amount of reasonable normalization can account for the 
excess of radio sources appearing for log.Li.4GHz > 24.5 — 
this would require a factor of at least 10 3 difference in 
the cluster and field LF normalizations, and the observed 
factor is ~ 180. 

5.3. Faint Cluster Sources 

The disparity between the faint-end slope of the IR clus- 
ter galaxies and the YRC field galaxies may indicate that 
(1) we do not have a good constraint on the sampled vol- 
ume for lower luminosity sources or (2) the cluster LF is 
truly flatter than the field LF at low luminosity. Explana- 
tion (1) can be illustrated by the observation that fainter 
sources in AWM7 are not detected as compared with A262 
and A426, which just encompass AWM7 in rcdshift space 
(see Figure 1). The statistics of AWM7 sources are poor, 
and the sampled volume may be much smaller than the 
one taken here. Statement (2) might be true if the num- 
ber counts of faint sources in clusters is lower than in the 
field due to the fainter and less massive sources accreting 
or merging with larger galaxies resulting in a top-heavy 
galaxy mass function. 

The IR LF alone indicates a flattening of the faint-end 
slope when compared with the field. This flattening is 
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Fig. 8. — The 1.4 GHz radio luminosity function for the 182 
sources in the composite cluster is shown by the open circles, and 
that for the late-type cluster galaxies is shown by the filled trian- 
gles. The solid curve denotes the best-fit field Schcchtcr LF for the 
radio data of the YRC field sample, normalized to the value of the 
cluster LF for the bin containing the greatest number of sources. 
We find the radio LF to more closely match that of the field sample. 
Three sources with logLi.4GHz > 24.5 are not shown here, and they 
contribute to the systematic excess of sources above the field LF for 
high luminosity objects, most of which is due to the radio luminous 
early-type galaxies. 



unlikely due to NED selection since our flux cutoffs are 
low enough to include normal field galaxies where optical 
and IR selection are identical (e.g., Condon et al. 2002, 
YRC). The effect responsible for the flattening of the faint- 
end slope (i.e., selection effects, top-heavy galaxy mass 
function, etc.) of the IR LF should also effect the radio 
LF faint-end slope in the same manner, when we take into 
account the radio- luminous Virgo sources with logLgo^m 
between 7 and 8. We do not currently have the means to 
ascertain the reason for such an effect, but merely suggest 
it for future study. 

6. DISCUSSION 

6.1. Density-Morphology Relation 

The density distribution computed from deprojected 
late-type Ring galaxies shows no trend with radius and 
is scattered around a mean of ~ 8 Mpc~ 3 . The computed 
density of late-type galaxies within the Core is only 2.0 
Mpc~ 3 , <~ 15 times under-dense than the composite clus- 
ter, and a factor of 14 smaller than the density of early- 
type galaxies in the Core region (28.4 Mpc -3 ). The over- 
all surface density of ellipticals in the Ring region is 2.86 
Mpc -2 , compared with a surface density of 19.10 Mpc~ 2 
in the Core region. These observations indicate that early- 
type galaxies are more strongly clustered in the Core than 
late-type galaxies. The bias of early-type galaxies and the 
analytic fit to the composite cluster density (§3) indicate 
that the high detection rate of ellipticals is not due to 
sample selection effects, but represents a true increase in 
density. 

The density-morphology relation is well established for 
optically-selected cluster galaxies, but ours is a highly bi- 
ased subset based on radio and FIR emission. A valid 
point of inquiry therefore is whether the increase in early- 
type galaxy density in the Core region is due to an under- 
lying change in the density-morphology relation for radio- 
and FIR-selected cluster galaxies, or if it simply reflects the 
existing density-morphology relation for all cluster galax- 
ies and a random, uniform appearance of radio sources 
among them. To answer this, we note that the exacting 
selection criteria used in this study makes it likely that 
the radio and FIR ellipticals in our sample are a subset 
of the optically-identified ellipticals. The total fraction 
of detected early-types (cD, elliptical, and SO) in both 
the flux- and luminosity- limited samples is ~ 30%, some- 
what lower than the fraction typically found in "spiral- 
rich" clusters (~ 50%). Consequently, the detection rate 
of early-type galaxies in our survey is consistent with a 
random population of radio sources tracing the underly- 
ing density-morphology relation for optical cluster galaxy 
samples. The distribution of late-type galaxies can be un- 
derstood in a similar manner to that of the early-type 
galaxies. Optically-selected samples show a late-type frac- 
tion that decreases toward the central regions of clusters. 
The same correlation is found in our sample (see above). 
Our flux cutoff is low enough to include galaxies with star 
formation rates similar to those found in the field, and 
the 60^im selection is identical to the optical selection of 
late-type galaxies (Condon et al. 2002, YRC). 

6.2. Mass Segregation in Clusters 
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Mass segregation among cluster galaxies may induce the 
observed correlations between radio and FIR luminosity 
and clustercentric radius as more massive galaxies may be 
more efficient IR and radio emitters. To investigate this 
possibility, if -band fluxes for all candidate sources were 
obtained from the Two Micron All Sky Survey (2MASS), 
using if-band luminosity as a proxy for mass. The if -band 
luminosity to mass ratio is thought to be more universal 
than optical ratios, subject to lower extinction effects, and 
linked to older stars. The 2MASS survey data include 
accurate magnitudes for sources less than 101" in extent. 

The cluster A1060, a few sources in AWM7 and Virgo, 
and almost all sources in Coma were observed by 2MASS, 
but have either not been cataloged, or, as in the case of 
Virgo, have not been mosaiced to obtain accurate 2.2/im 
fluxes. For example, spatially integrated 2.2/im fluxes are 
not reported for unmosaiced Virgo cluster galaxies greater 
than 101" in extent. These are some of the most massive 
galaxies in our sample. To assess possible mass segregation 
among Virgo galaxies, we have used the data obtained by 
Boselli et al. (1997), which includes if' magnitudes mea- 
sured from Infrared Space Observatory (ISO) images of 
102 Virgo cluster galaxies. The data for Virgo are an- 
alyzed seperately from the remaining sample to preclude 
any systematic bias owing to differences in wavelength cov- 
erage and response between the 2MASS if-band and the 
ISO if '-band filters. 

For the 54 sources with reported 2MASS if-band mag- 
nitudes, the if-band absolute luminosity has been com- 
puted and plotted against projected radial distance (Fig- 
ure 9a). A formal fit shows only a very weak correlation 
with a = 1.11 dex: a ~ 45% probability of a positive cor- 
relation between M(K) and r. Separate \ 2 analyses for 
early- and late-type galaxies find the correlation is even 
weaker: the corresponding errors in the slope of the corre- 
lations allow for constant M(K). The slight correlation for 
the composite cluster may be due to the central cD galax- 
ies, as well as the Messier and Seyfert galaxies in Virgo 
with an extent of < 101" in diameter. The formal fit 
was redone excluding the two cD galaxies (NGC 708 and 
NGC 1275) for which if-band data were available. Some 
galaxies with very small cluster offsets, small velocity dif- 
ferences, or unusually low q value were left in the sample if 
their optical morphology indicated a non-cD galaxy. The 
resulting correlation is even weaker (Figure 9a), and the 
slope is statistically consistent with zero (0.24 ± 0.39 with 
a = 1.07 dex). Eleven galaxies in our Virgo sample were 
imaged in the Boselli et al. (1997) study, including half 
of our detected Messier sources. Their absolute if' -band 
magnitudes are plotted in Figure 9b as a function of pro- 
jected radial distance. As with the previous analysis, we 
find only a ~ 45% probability of a positive correlation be- 
tween M(K) and r, and the statistics are consistent with 
a flat slope. 

Luminosity segregation in clusters has been studied ex- 
tensively, leading to mass estimates by invoking some M /L 
ratio. The mass per galaxy can be obtained by normal- 
izing the luminosity density by the number density. As 
an example, we consider the Virgo cluster. The luminos- 
ity density profile has an identical slope to the number 
density distribution in Virgo, except in the very central 
(r < 60') region of the cluster, where the luminosity den- 
sity is steeper than the number density (Schindler et al. 
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Fig. 9. — (a) Absolute A"-band magnitude as a function of radius 
for the sample of 54 galaxies with reported 2MASS if-band fluxes. 
Filled circles are early-type galaxies and open squares are late-type 
galaxies. The solid and dotted lines indicate the best-fit correlations 
for the sample with and without cD galaxies, respectively, (b) Same 
as (a), but for the 11 Virgo cluster galaxies in the Boselli ct al. 
(1997) near-IR imaging survey. In both samples, only a very weak 
positive correlation is found, indicating that massive galaxies in the 
composite cluster do not preferentially populate cluster cores. 



1999). Neglecting these central bright galaxies (most of 
which are Messier and Seyfert galaxies), the inferred av- 
erage mass per galaxy is constant with radius. Again, 
this conclusion assumes some constant M/L ratio based 
on the average M/L ratio of early- and late- type galaxies 
in clusters. As another example, there is marginal evi- 
dence for segregation in the inner 20% of the core radius 
of Coma, but no large-scale correlation is found (Zhao et 
al. 1991; Bahcall 1973 and references therein). Inclusion 
of Coma galaxies in Figure 9 should therefore not signifi- 
cantly affect our conclusion based on the clusters for which 
if-band data were available. The absence of mass segre- 
gation for the time-averaged evolution of the composite 
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cluster is important for our study since it indicates that 
radial variations in radio and FIR emission are not a result 
of gas-poor massive ellipticals dynamically sinking into the 
Core region. 

6.3. Discerning AGN from the Normal Population 

Determining the power source of cluster galaxies is cru- 
cial in order to correctly interpret the radio and FIR corre- 
lations discussed in previous sections. Our sample is small 
enough that a case-by-case analysis of the multiwavelcngth 
data for each galaxy can be carried out to determine which 
show evidence of AGN. We have searched the literature for 
each of the radio-excess sources, and we assign an AGN 
classification based on evidence of hard X-ray emission, 
radio jets, or broad optical emission lines consistent with 
an active nucleus. This resulted in 5 AGN classifications. 
Miller & Owen (2001a) obtained high- resolution data for 
two of our 4 remaining radio-excess sources (CGCG097- 
073 and KUG1258+287) and were classified as star form- 
ing galaxies based on their nuclear optical spectra. 

The validity of using radio-excess as an indicator of AGN 
has been explored in previous studies (e.g., YRC; Yun et 
al. 2000; Condon & Broderick 1988). The exact AGN 
fraction among radio-excess sources in the Miller & Owen 
(2001a) sample is unclear as (1) they were unable to obtain 
spectroscopy for all of their sources, (2) their spectroscopy 
cannot distinguish a weak AGN system from a compact 
starburst dominated emission-line spectrum, and (3) they 
include galaxies in the AGN sample which show evidence 
of an old stellar population from the absorption-line spec- 
trum (Miller & Owen 2001a). They test the effects of 
(2) by obtaining high-resolution (~ 0.25") radio images 
of cluster galaxies in A1367 and A1656, and find that of 
the 5 of 20 detected galaxies, at least 4 have a compact 
nucleus. Of the remaining 15, 13 are known star form- 
ing galaxies whose nebular emission may overwhelm the 
emission associated with an AGN. 

We have independently analyzed the data of Miller & 
Owen (2001a) and find that of the ~ 70 spectroscopically- 
confirmed star forming galaxies in their radio-selected 
sample, only 6 show radio-excess, whereas ~ 23/77 = 30% 
of their spectroscopically-confirmcd AGNs show radio ex- 
cess. The overall frequency of AGN sources among their 
radio-excess galaxies is therefore ~ 80%. It is likely that 
the actual frequency of AGNs among the radio-excess 
galaxies is higher, since few optical spectra were collected 
for the most radio-luminous sources (q < 1.00). Miller 
& Owen (2001b) detect 4 FIR luminous (logipm > 11.0) 
galaxies. Of these, one is an IR-excess galaxy with q = 3.21 
and the rest have nominal q values. None of the IR-excess 
sources in our sample show evidence for AGN activity. The 
13 radio-excess sources in our sample include 4 AGN cD 
galaxies and 5 non-cD AGN sources. Hence, of the radio- 
excess sources in our sample, ~ 70% are AGN based on 
the presence of radio jets, X-ray emission, and/or optical 
emission lines. This fraction is a lower limit since some of 
the less well-studied radio-excess objects could also host 
an AGN. In summary, it appears that a significant frac- 
tion of radio-excess objects are associated with luminous 
AGN. 

A segregation between IR luminous, normal, and radio- 
excess objects can be seen by examining the IRAS FIR- 
vs.-mid-IR color-color diagram (see YRC for a discussion) . 



The comparison of the environment in radio-excess cluster 
objects to "normal" cluster galaxies can be investigated in 
the same way. In Figure 10, the large open and filled grey 
circles represent early- and late-type cluster galaxies with 
q > 1.64, respectively. The radio-excess sources are shown 
with large open and filled squares for late- and early-type 
galaxies, respectively. For comparison, the distribution 
for the IR-selected YRC field galaxies within 94 Mpc is 
indicated by the small black points. 
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Fig. 10. — IRAS FIR versus mid-IR color-color diagram for those 
sources in the luminosity-limited sample with available IRAS 25/im 
measurements. FIR color has been converted to temperature using 
j8 = 1. The large open and filled grey circles represent early- and 
late-type galaxies, respectively. The radio-excess sources are shown 
with large open and filled squares for early- and late-type galaxies, 
respectively. The small points are for the YRC field galaxies with 
distance less than 94 Mpc. The hashed vertical line denotes the 
proposed line indicating the lower boundary for IR Seyfert galaxies 
(de Grijp ct al. 1985), and radio-excess sources are found to the right 
of this line. This suggests that the mid-IR excess which characterizes 
these sources is attributable to dust heating in the circumnuclear 
environment surrounding an AGN. 

The aggregate sample of cluster galaxies in the luminosity- 
limited sample have (S^opm/ Sioo^m) ~ 0.42±0.17, and the 
appropriately selected field galaxies have (Seo^m/Sioo/jm) ~ 
0.49 ± 0.07. For comparison, the mean FIR color in the 
cluster sample of Bicay & Giovanelli (1987) is (S , 6 o M m/'S , ioopm) 
0.38, slightly lower, but still consistent with our result. 
The slight difference may be due to our selection of clus- 
ter galaxies based on FIR emission, biasing the sample 
to galaxies with warmer FIR color than seen in optically- 
selected galaxies. 

All of the radio-excess sources which have detected 
25/im emission appear to the right of the proposed line 
indicating the lower boundary for IR Seyfert galaxies ac- 
cording to their S^m/SWiri colors (de Grijp et al. 1985). 
Sources detected at 25yum will have higher <S , 25 A mi/5Wm 
ratios in general than those not detected at 25^m. The 
absence of radio-excess sources to the left of the proposed 
indicator for IR Seyfert galaxies is not significant by itself. 
However, when taken in conjunction with the observation 
that most radio-excess galaxies in the field have warmer 
mid-IR colors than the normal population by ~ 0.5 dex 
(YRC), this suggests that radio-excess sources in the clus- 
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ter also, on average, have warmer mid-IR colors, consistent 
with a model in which dust is heated in the circumnuclcar 
region by an AGN. The early- type radio-excess objects 
(most of which are AGN) populate the same region in 
color-color space as the radio-excess field objects appear- 
ing to the right of the IR Seyfert line (YRC). In summary, 
radio-excess sources show mid-IR colors consistent with 
AGN, and this is supported by the case-by-case analysis. 
The segregation of q values indicates that most of these 
radio-excess sources are in the Core region. 

While radio-excess sources generally show mid-IR ex- 
cesses, the opposite is not true. Namely, there are many 
late-type galaxies that do not have a radio-excess but 
which do show mid-IR-excess. This is not surprising since 
cither a young or compact starburst or radio-quiet AGN 
could result in the observed mid-IR color (e.g., Efstathiou, 
Rowan- Robinson, & Siebenmorgen 2000). We can place an 
upper limit to the fraction of mid-IR excess cluster sources 
with q > 1.64 of 51%, with the field value being < 10%. A 
direct comparison of the two fractions is risky as there may 
be a number of selection effects we have not considered. 
If, in fact, a larger fraction of normal cluster galaxies ex- 
hibit mid-IR excess when compared with the field, it may 
indicate that many of the disk galaxies host young, dust- 
enshrouded starbursts. If so, this provides some marginal 
evidence for harassment-induced star formation activity 
(see §6.6). 

6.4. FIR Emission Mechanisms 
6.4.1. X-Ray Dust Heating 

According to §6.2, the composite cluster shows no evi- 
dence of mass segregation based on the galaxies' if-band 
luminosity. If only the cirrus component of dust is stripped 
from a galaxy, then this will have an effect on the observed 
emission at lOO^im, but the net impact on Lfir should be 
small. This argument is difficult to quantify since it re- 
quires a knowledge of the dust distribution within galax- 
ies. However, dust is closely associated with molecular gas 
which can be accurately traced by millimeter CO tran- 
sitions. Young et al. (1995) have shown that more than 
90% of the CO emission from late- type galaxies comes from 
the inner 2/3 of the optical disk, and the same may also 
be true of FIR emission. Consequently, the mass-to-light 
ratio as traced by M(K)/L 60Aim should be insensitive to 
the effects of dust stripping. The presence of significant 
collisional heating should result in the more centrally lo- 
cated galaxies having lower M(K)/L6 UM m values. Figure 11 
shows no correlation between M(K)/L 60pm and radius in- 
dicating that collisional heating due to electrons is not a 
significant source of FIR emission. 

The effects of collisional heating of dust by ICM elec- 
trons can be estimated by calculating the collisional heat- 
ing rate. Following the analysis of Goudfrooij & de Jong 
(1995), the heating rate is given by 



H e (r) = n c (r) 



8kT c 



1/2 



-kT r 



ira 



(9) 



where T e is the electron temperature, a is the dust grain 
radius, and n c (r) is the radial electron density profile, as- 
sumed to follow a f3 profile similar to equation (5). The 
dust temperature is obtained by equating the collisional 
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Fig. 11. — Ratio between _R"-band absolute magnitude and log of 
60/im luminosity versus clustercentric distance for the luminosity- 
limited sample. Solid circles are early-type galaxies and open 
squares are late-type galaxies. We find no significant correlation 
between the two, indicating that electron-dust collisions have a rel- 
atively insignificant effect on the dust heating rate for most galaxies 
in the cluster environment (see text). 



heating rate to the cooling rate due to IR emission given 
by de Jong et al. (1990), yielding 



n e (r)(8fcT e Am e ) 1 / 2 (3fcr c /2) 
6.0 x 10- 8 K o (4ap d /3) 



1/5 



(10) 



where the opacity n — 2.4 x 10 3 (Hildcbrand 1983), pd 
is the grain density, and a dust emissivity proportional 
to A -1 has been assumed. Assuming a typical value of 
T ~ 10 7 K, the temperature-radius curves for 6 of the 7 
clusters (excluding Virgo) are shown in Figure 12a, for the 
values of r c , (3, and n o (0) derived by Mohr, Mathiescn, & 
Evrard (1999) from the X-ray surface brightness profiles. 
The curves shown are those corresponding to a = 10~ 5 
cm and pd — 3 g cm~ 3 . The temperature-radius relation 
is sensitive to the product of a and pd, and Figure 12b 
shows the curves corresponding to different grain densities 
assuming the /3-model parameters derived in §3 for the 
composite cluster. 

The effects of X-ray dust heating do appear to be signif- 
icant in cD galaxies containing large amounts of hot X-ray 
gas associated with the ICM . In this case, the electron 
density is high enough that dust heating becomes impor- 
tant (e.g., Bregman, McNamara, & O'Conncll 1990; Lester 
et al. 1995) and the emission is uncharacteristic of the lo- 
calized FIR emission associated with active star formation. 
This can be seen in Figures 12a and 12b, where the dust 
temperature at the center of the cluster is ~ 40K, 21% 
higher than the average temperature derived for cluster 
galaxies using the FIR color alone. 

Figure 11 indicates that collisional heating of dust grains 
in the sample galaxies has a negligible effect on the ob- 
served FIR emission, allowing us to place a loose constraint 
on the dust column density. For example, the mean FIR 
color temperature for the cluster galaxies is ~ 33K (see 
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Fig. 12. — (a) Dust temperature versus radius for a collisionally- 
heated grain population of size a = 10 — 5 cm and grain density 
pd = 3 g cm -3 for 6 of the 7 sample clusters, using the /3-model 
parameters derived from the X-ray surface brightness profile (Mohr, 
Mathicsen, & Evrard 1999). (b) Dust temperature versus radius 
for a collisionally-hcatcd grain population of size a = 10~ 5 cm for 
the composite cluster, adopting the /3-model parameters from §3. 
The solid and hashed curves correspond to p d = 1 and 3 g cm -3 , 
respectively. The solid vertical lines in Figure 12a and Figure 12b 
denote the boundary between the Core and Ring regions. 

Figure 10), indicating a/?d > 2.7 x 10~ 5 g cm~ 2 . For the 
typical values of a and in Figure 12b, the derived dust 
temperatures are not sufficient to account for the observed 
FIR color temperature, particularly in the Ring region, 
where Td falls below 30 K. If we assume a limiting value 
of apd = 2.7 x 10~ 5 g cm~ 2 , then the collisional heating 
rate is <J 50% of the optical heating rate at r = 0.5 Mpc 
and <; 6% of the optical heating rate at r = 1.5 Mpc. The 
lower limit of apd given above indicates that large dust 
grains are primarily responsible for the FIR cmissivity. 



6.4.2. Gas and Dust Stripping 

Spiral galaxies in rich clusters show a marked HI gas 
deficiency when compared with isolated field spiral galax- 
ies (Giovanelli & Haynes 1985; Cayette et al. 1990). The 
mechanisms used to explain stripping include ram pres- 
sure (Gunn & Gott 1972), turbulent viscosity (Nulscn 
1982), tides during galaxy-galaxy interaction (e.g., Rich- 
stone 1975; Merrit 1983), or thermal evaporation (Cowie 
& Songaila 1977). More recent studies have suggested that 
the inner disk of molecular gas remains intact as it is more 
strongly bound to the galaxy (Kenney & Young 1986; Cay- 
atte et al. 1994; Perea et al. 1997; Bicay & Giovanelli 
1987). Contursi et al. (2001) have found that HI deficient 
galaxies in Coma and A1367 generally have comparable 
dust masses to HI- normal galaxies. 

The FIR emission for cluster galaxies in our sample is in- 
dependent of radial distance in agreement with that found 
in other cluster samples (e.g., Rengarajan, Karnik, & Iyen- 
gar 1997). The apparent flatness of the FIR luminosity- 
radius correlation does not address directly the question 
of the gas and dust content of the galaxies since a sys- 
tem with globally low gas content could have localized star 
formation producing the same levels of FIR emission as a 
gas-rich system with modest levels of global star forma- 
tion. While higher resolution data can better address this 
question, optical (e.g, Ha), mm (e.g., CO), and radio (e.g., 
HI 21 cm) tracers indicate that galaxies have lower cold gas 
content toward the central regions of clusters. This trend, 
if it is exists in our sample, does not affect the observed 
FIR emission. 

The flat FIR luminosity-radius correlation seen in Fig- 
ure 6 suggests that changes in the environment within 1.5 
Mpc do not have a significant impact on the FIR emission. 
The cluster-to-cluster variation in the FIR luminosity- 
radius correlation is shown in Figure 13 for 6 of the 7 clus- 
ters. A262 and AWM7 have few detected sources above 
our luminosity cutoff, and any trend will be unclear given 
the few data points for these two clusters. For the re- 
maining clusters, the FIR luminosity is relatively flat with 
radius. 

The distinct lack of detected radio and FIR sources in 
the Virgo cluster core could be a result of severe FIR sup- 
pression (Figure 14a). This reduction in FIR luminosity 
is correlated with HI deficiency in the Virgo galaxies stud- 
ied by Doyon & Joseph (1989) — many of which fall within 
our Core region — and who quote a factor of 3 decrease in 
star formation as a result. The recent Ha imaging study of 
Virgo spiral galaxies by Koopmann & Kenney (2002) indi- 
cates that almost all show evidence for truncated star for- 
mation. The FIR emission for Ring Virgo galaxies appears 
to weakly decrease with radius (Figure 14a), and may be 
expected if the diffuse cold dust component is stripped in 
the higher density regions of the cluster, leaving only the 
warm component to dominate the FIR luminosity. To test 
for this effect, we compute S , 60pm/<S'ioopm (FIR color) ver- 
sus radius for the flux-limited sample of Virgo late-type 
galaxies (Figure 14b). Stripping of cold dust should yield 
warmer FIR colors toward the cluster center, but this is 
not seen in Figure 14b. We conclude that cold dust strip- 
ping, if it exists in the Virgo cluster, does not affect the 
FIR emission. 
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Fig. 13. — IR luminosity versus projected radius for 6 of the 7 clusters in the flux-limited sample. The horizontal line in each plot 
denotes the completeness limit. Filled circles and open squares represent early- and late-type galaxies, respectively. We find no significant 
clustcr-to-clustcr variation in the IR luminosity as a function of radius. 



6.5. Radio Enhancement 

The q value is not affected by any increase or decrease 
in the star formation rate of a galaxy as the radio and FIR 
emission are directly proportional to this rate. Our obser- 
vation of a g-radius correlation and lack of FIR luminosity- 
radius correlation (§4.3) indicate that changes in star for- 
mation rates are not responsible for the increase in radio 
luminosity of sources toward the center of the composite 
cluster. Tidally-induced star formation is therefore not a 
dominant mechanism in increasing the radio emissivity of 
sources toward the center clusters. 

While there are relatively few mechanisms that can only 
affect the observed FIR emission, radio enhancement with 
nominal FIR emissivity can occur in a number of ways, 
including ram-pressure and thermal magnctospheric com- 
pression, £ md elevated AGN activity. Each of these possi- 
bilities is discussed below. 

6.5.1. Ram Pressure Magneto spheric Compression 

Ram pressure magnetospheric compression can occur if 
a galaxy has sufficient velocity through the ICM causing 
the gas on the leading edge of the galaxy to be compressed. 
The magnetic field coupled to the gas will also be com- 
pressed, resulting in increased synchrotron emission. Al- 
ternatively, thermal compression occurs when the pressure 
of the ICM itself is enough to compress the magnetic field 
to the point where it may affect the observed synchrotron 
emission (see below). 

The effects of ram pressure are assessed by computing 
the line-of-sight velocity deviation of each galaxy from the 
cluster velocity normalized by the cluster velocity disper- 



sion: 

\v galaxy ^cluster | (11) 

0~ cluster 

referred to as the "relative velocity" of each galaxy (see 
Miller & Owen [2001a], eq. [3]). We have obtained the 
cluster velocity dispersion data for the Abell sources from 
Struble & Rood (1999). The velocity dispersion of AWM7 
and Virgo were obtained from Koranyi et al. (1998) and 
Fadda et al. (1996), respectively. The velocity disper- 
sion of AWM7 shows a weak radial profile, with dcx ~ 150 
km s _1 lower in the inner 0.13 Mpc. Since our sample ex- 
tends to 1.5 Mpc, we adopt the dispersion velocity found 
for the outskirts of the cluster (a ~ 680 km s _1 ). A sum- 
mary of the adopted cluster velocity dispersions is given 
in Table 1. 

Figure 15 shows the distribution of q values as a function 
of relative velocity for the cluster galaxies, with separate 
emphasis on early and late-type galaxies. The deficit of 
galaxies in the lower right-hand corner of Figure 15 is due 
to projection effects, since some of those galaxies with a 
large cluster velocity may have small line-of-sight velocity 
if the direction of motion is aligned perpendicular to the 
line-of-sight. 

A significant contribution of ram pressure magneto- 
spheric compression should manifest itself as an anti- 
correlation between the q distribution and relative veloc- 
ity. For early-type galaxies, we find a slope of 0.56 ± 0.28 
(including radio AGNs). Excluding all early- type galax- 
ies with radio-excess, the slope becomes zero. One might 
argue that excluding all radio-excess sources is not ap- 
propriate since some of the radio-excess in these galaxies 
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Fig. 14. — (a) IR luminosity versus projected radius for the flux- 
limited sample of Virgo galaxies. Filled circles and open squares 
indicate early- and late- type galaxies, respectively, and the solid 
horizontal line denotes our luminosity limit, (b) IRAS FIR color as 
a function of projected radius for the flux-limited sample of Virgo 
galaxies. Filled circles and open squares indicate early- and late- 
type galaxies respectively. There is no correlation between FIR 
color and radius for both early- and latc-typc galaxies. Also note 
the paucity of IR sources in Virgo for radii < 0.5 Mpc from the 
cluster center. 



may be due to ram pressure magnetospheric compression, 
particularly in light of the efficiency for ram pressure to 
increase synchrotron emission due to the v 2 dependence 
of the magnetic field. However, most of these sources do 
show evidence of powerful AGN activity and most are mas- 
sive ellipticals with rather low relative velocity (Figure 15). 
Projection effects may also produce this effect, but it is 
unlikely that all radio-excess sources have velocities per- 
pendicular to the line-of-sight. Therefore, the removal of 
these radio-excess sources is warranted given the appar- 
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Fig. 15. — q parameter as a function of relative (projected) ve- 
locity, as computed from equation (11), for the luminosity-limited 
sample. Filled circles and open squares represent early- and late- 
type galaxies, respectively. We find only a very weak correlation 
between q and 5v. The correlation, if it exists, may be masked by a 
number of effects (see text). 



cntly distinct source for radio-emission. Hence, there is 
no obvious trend indicating that increased velocity among 
elliptical galaxies results in increased radio emission. 

One caveat is as follows. Low velocity ellipticals in 
dense environments will have a longer interaction time 
with neighboring galaxies. The result may be increased 
gas accumulation in the major component of the interac- 
tion, leading to an increase in the intensity and/or dura- 
tion of star formation activity and/or result in elevated 
AGN activity. Our data suggest that tidally-induced star 
formation is not a major contributor to the enhancement 
of radio emission in most non-cD elliptical galaxies. As- 
suming this to be the case, it is unlikely that a longer 
interaction time is important in increasing q values for 
low-velocity galaxies in the Core region. The low-density 
environment characteristic of late-type galaxies (cf. §3) in 
the Ring argues against the possibility of increased radio 
emission due to longer interaction time: the velocity dis- 
persion generally flattens out for large radii, whereas the 
density falls off rather rapidly as a power law (cf. eq. [5]). 

The late-type galaxies exhibit a weak correlation be- 
tween q and Sv, similar to the findings of Miller & Owen 
(2001a). Other studies have indicated little correlation 
between projected velocity and projected radius. As an 
example, there is great scatter between velocity and ra- 
dius in the Coma cluster with no obvious relation between 
the two (Biviano et al. 1995). The statistics of our sample 
are such that we cannot make a secure determination for 
individual clusters, and the lack of a strong q-Sv correla- 
tion does not rule out the possibility that some of the radio 
emission from the cluster galaxies is due to ram-pressure 
enhancement of the magnetic field. 

There are clear cases where ram pressure enhancement 
has a noticeable effect on the synchrotron power (e.g., 
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Cayette et al. 1990). However, the effect, even if present 
on a large scale, may be difficult to observe. In particu- 
lar, Miller & Owen (2001a) suggest several possibilities for 
non-gaussian motions in clusters, including bulk flows in 
clusters and cluster cooling flows. Projection effects may 
also weaken such a correlation if it exists. These effects, 
compounded with a lack of knowledge of the spatial coor- 
dinates of a galaxy with respect to the core of the cluster, 
all serve to muddle any correlation between q and Sv. 

6.5.2. Thermal Compression 

The observed radio enhancement among star forming 
late-type galaxies could result from magnetic field com- 
pression due to the thermal pressure of the ICM. Miller & 
Owen (2001a) find that thermal pressure can account for 
most, if not all, of the pressure required to produce the ob- 
served radio enhancement in their sample. The effects of 
thermal compression can be quantified by computing the 
increasing in synchrotron power due to the compression of 
the magnetic field as a result of external pressure. Follow- 
ing the analysis of Miller & Owen (2001a), the gas density 
distribution will follow the (5 profile of equation (5), and 
the thermal pressure is 



Pth 



pkT 
pm R ' 



(12) 



where T is the temperature of the gas and p, = 1 for a gas 
of pure hydrogen. The synchrotron power is directly pro- 
portional to the pressure on the gas which, in the cluster, 
can be thought of as the sum of two pressures: 



Pot = P{B ) + Pth = 



IP 
8tt 



(13) 



Here, P(B Q ) refers to the pressure due to the ambient mag- 
netic field of the galaxy if it were not located in the cluster, 
Pth refers to the thermal pressure provided by the ICM, 
and B refers to the value of the magnetic field that would 
produce such a total pressure. The ratio of the synchrotron 
luminosity in the cluster (L c ) to that of the same galaxy 
in the field (Lf) can then be written as 



Lc 
U 



P(B ) + P th 
P(Bo) 



1 



pkT 



^m H P(P ) 



(14) 



The dashed curve in Figure 6b represents the expected 
value of L c as a function of radius, where we have made 
the following assumptions. First, similar to Miller & Owen 
(2001a), we have assumed the /3-model parameters for 
Coma as deduced by Mohr et al. (1999), as it is most rep- 
resentative of relaxed clusters (r c — 0.257 Mpc, — 0.705, 
p c = 7.4 x 10~ 27 g cm~ 3 , and kT = 8.21 keV). Secondly, 
we have assumed that a typical field spiral galaxy has a 
magnetic field of ~ bpG. For these values, the thermal 
pressure is ~ 1.35 x 10~ 12 dyn cm~ 2 at r = 1.5 Mpc from 
the cluster center, approximately 1.35 times that of the 
magnetic pressure associated with a 5/xG field. Equation 
(14) can be normalized by assuming that the mean radio 
luminosity at r = 3.0 Mpc is L 1ACmz = 10 20 61 W Hz" 1 , 
according to the best-fit correlation (solid line) shown in 
Figure 6b, and that this value is comparable to the lumi- 
nosity of field galaxies. 



For this analysis, the model of radio luminosity expected 
due to thermal compression is in general agreement with 
the observed data. However, the dashed curve in Figure 6 
should only be considered as representative of the effects 
that thermal pressure has on the synchrotron power. In 
particular, the normalization of the curve is subject to a 
number of factors, including the assumed B field for spiral 
galaxies and the assumed (3 model parameters. The ex- 
trapolation of the best-fit curve between radio luminosity 
and radius for the cluster sample out to 3 Mpc may also 
be inaccurate as we have no data that samples the field 
environment around the clusters. We also note that the 
intrinsic scatter in the data makes it difficult to determine 
if the underlying correlation is truly nonlinear, as expected 
for thermal pressure effects, or if it is linear, as it appears 
in Figure 6b. Note that the density used in equation (14) 
is a deprojected quantity. Therefore, the dashed curve in 
Figure 6b increasingly overestimates the radio luminosity 
for decreasing radius. Correcting for this effect will bring 
the dashed curve downward toward the solid line for small 
radii (e.g., within the Core region). In summary, we find 
that thermal compression could explain the relative change 
in the radio luminosities of cluster galaxies located within 
1.5 Mpc of cluster centers. 

6.6. Delayed Harassment and Frequency of Radio-Excess 

Sources 

A recent model of galaxy-galaxy interaction has been 
proposed by Lake, Katz, & Moore (1998). They carry 
out numerical simulations of disk galaxies in rich clus- 
ter environments and show that galaxy "harassment" by 
many galaxy-galaxy interactions can drive up to 90% of 
the gas into the inner 500 pc. This gas infall can lead 
to radio enhancement by (1) fueling AGNs (see §6.5.3) 
and/or quasars in low-luminosity hosts (Lake, Katz, & 
Moore 1998), or (2) nuclear starburst activity (e.g., Fujita 
1998). As discussed previously, we can rule out (2) being 
the sole mechanism in producing the correlations discussed 
in §4.3 since it only affects the star formation rate. 

While harassment is dominant in only the densest clus- 
ter environments (e.g., cluster cores), the effects of harass- 
ment may also influence the Ring population as they may 
not fully materialize until some time has passed. If AGN 
activation is immediate due to gravitational, hydrody- 
namic, or magnetodynamic perturbation, then we should 
expect sources with lower q values to have the largest rel- 
ative velocity. Figure 15 suggests the opposite is true, but 
testing this would require disentangling projection effects 
(see §6.5.1). 

We note that the Core crossing time is ~ 1.7 Gyr for a 
galaxy with a relative cluster velocity of ~ 600 km sec -1 
(similar to the velocity dispersion), much longer than the 
disk crossing time of about 100 Myr for a starburst. Con- 
sequently, simple delayed harassment is an unlikely cause 
for the observed correlation between q and r, as the in- 
duced activity should dissipate by the time galaxies enter 
the Ring region. A more complicated model of harassment 
(e.g., such as one that includes the presence of cluster sub- 
structure) may be needed to explain both the correlation 
between q and r and the wide-spread behavior of the mid- 
IR sources. 

The radio-excess fraction among cluster galaxies is 11 
times higher than that among IR-selected field galaxies. 
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The sources in the spectroscopically confirmed AGN sam- 
ple of Miller & Owen (2001a) show a correlation between q 
and r similar to that of the star forming galaxies, suggest- 
ing that some of the radio emission in AGN sources is due 
to normal star formation in the disk of the galaxy (Roy et 
al. 1998; Baum et al. 1993). It is therefore plausible that 
the enhancement in radio emission due to thermal pres- 
sure is sufficient to increase the radio luminosity of AGN 
hosts such that they fall above the threshold to be con- 
sidered a radio-excess object. In particular, the effects of 
thermal pressure result in a mean radio luminosity that 
is 4 times greater in the Core than in the Ring region, 
implying a reduction in q value of ~ 0.6 dex. The 6 radio- 
excess sources with q > 1 may not have been observed as 
radio-excess sources if they were in the field, implying an 
observed radio-excess fraction of 3/114 ~ 2.6% (excluding 
cD galaxies), which is still larger than, but consistent with, 
the radio-excess fraction found by YRC (~ 1.3%). Note 
that despite the possible effects of radio enhancement via 
thermal compression in our sample, the large majority of 
the observed radio-excess sources (70%) are AGNs. 

7. CONCLUSIONS 

We have studied the effects of the cluster environment 
on the star formation and AGN activity in cluster galaxies. 
To accomplish this, radio and FIR data were compiled for 
a luminosity-limited sample of 114 cluster galaxies in the 
most massive, nearby (< 100 Mpc) X-ray clusters. The 
selection based on X-ray luminosity ensures that we are 
including only the most massive and dynamically-related 
clusters, and is a direct probe of the state of the intra- 
cluster gas. The seven clusters studied are A262, AWM7, 
A426, A1060, A1367, Virgo, and Coma. The main conclu- 
sions from this study are as follows: 

1. The radio-FIR correlation for the luminosity- limited 
sample shows almost three times the scatter in the correla- 
tion for field galaxies. An examination of the logarithmic 
FIR-radio flux density ratio (q) reveals that this scatter is 
mainly due to radio-excess (q < 1.64) sources, ~ 70% of 
which show evidence of AGN activity. Including these ex- 
treme sources, we find {q} = 2.07. The larger cluster sam- 
ple of Miller & Owen (2001a) has (q) = 1.95, lower than 
our value, and is likely due to their selection based strictly 
on radio emission. Similar field samples show (q) = 2.34 
(YRC). 

2. The radio-excess fraction in the Core region of the 
composite cluster (R < 0.5 Mpc; 21%) is larger than that 
found in the Ring region (0.5 < R < 1.5 Mpc; 7%), which 
is in turn larger than that found in the field (1%), indicat- 
ing a correlation with cluster position. We find a positive 
correlation between q and r (clustercentric distance), and 
this is a result of radio enhancement toward the center of 
clusters. The Miller & Owen (2001a) sample also indicates 
enhanced radio emission by a factor of 2 — 3 among core 
galaxies. 

3. We have computed the IR and radio cluster LFs. The 
excess of cluster sources above the field radio LF at high 
luminosity, as well as the excesses observed for the YRC 
sample and UGC sample all suggest that radio-luminous 
AGNs contribute more to the radio luminosity density in 
the cluster environment than they do in the field. 

4. The radio and morphological distributions indicate 
that radio-luminous and early-type galaxies, respectively, 



are clustered in the Core region. We find no evidence that 
the increase in density is due to a change in the underly- 
ing density-morphology relation due to our radio and FIR 
selection criteria. Rather, we are tracing only a subset 
of optically-identified early-type galaxies. Consequently, 
the increase in early-type density in the Core region re- 
flects the existing density-morphology relation for all clus- 
ter galaxies. 

5. As more massive galaxies may be more efficient ra- 
dio and IR emitters, we have investigated the possibility 
of mass segregation in our sample by compiling X-band 
magnitudes for all candidate sources. Based on the analy- 
sis of the time-averaged evolution of the composite cluster, 
we find that massive galaxies are not preferentially clus- 
tered in the Core region. It is therefore likely that the 
radio-luminous elliptical galaxies were formed in the Core 
region. 

6. An examination of the IRAS colors of the cluster pop- 
ulation indicates that radio-excess sources have warmer 
mid-IR colors than the normal cluster and field popula- 
tions, consistent with a model in which dust is heated in 
the circumnuclear region by an AGN. We have indepen- 
dently analyzed the Miller & Owen (2001a) radio-selected 
sample and find that the frequency of AGNs among their 
radio-excess sources is ~ 80%. 

7. ICM electron-dust collisions do not play a significant 
role in heating the dust in most cluster galaxies. The flat 
FIR luminosity-radius relation suggests that environmen- 
tal effects within the virial radii of clusters do not have a 
significant impact on the FIR emission, despite the factor 
of 10 decrease in density from the center of the cluster to 
1 virial radius (~ 1.5 Mpc). 

8. Among the mechanisms resulting in increased ra- 
dio emission, we cannot rule out a contribution from ram 
pressure enhancement of synchrotron emission. We find 
that thermal compression of a galaxy's ISM by the intr- 
acluster medium could also explain the relative enhance- 
ment in the radio luminosity of centrally-located galaxies 
compared with those at 1.5 Mpc. Thermal compression 
may increase the synchrotron emission from the disk of 
a galaxy and push radio-luminous AGNs into the radio- 
excess regime. We argue based on a timescalc argument 
that simple delayed harassment is an unlikely cause for 
the observed correlations, and a more complicated model 
including cluster substructure may be needed. 
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